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^ ■ Abstract 

\ In this paper, we investigate the stability of cognitive cooperative systems. We consider and propose 

three systems. The first system is the cooperative secondary transmitter-receiver system (CSTR), where, 
the secondary transmitter and the secondary receiver increase the spectrum availability for the secondary 
C/3 ' transmitter packets by relaying the unsuccessfully transmitted packets of the primary transmitter. We 

l— 'I provide inner and outer bounds for the CSTR system. The second system is the non-cooperative 

secondary receiver with no priority (NCNPQ) queueing to any of the queues. In the second system only 

> ■ 

the secondary transmitter helps the primary user to deliver its undelivered packets by certain access 

O ■ 

' probabilities assigned to each queue. The third system is the non-cooperative secondary receiver with 

' priority queueing (NCPQ), where, the secondary transmitter assigns priority queueing to the relaying 

packets. The maximum stable throughput regions of the considered systems are characterized. The 
results reveal that the maximum stable-throughput region of the NCPQ is a subset of the maximum 
j> ' stable-throughput of CSTR and NCNPQ. Also, the maximum stable-throughput of the NCNPQ is a 

' subset of the maximum stable-throughput of the CSTR system. 

Index Terms 

Cognitive radio, convex optimization, closure, stability analysis. 



I. INTRODUCTION 

The electromagnetic radio spectrum is a precious resource, the use of which by transmitters and receivers is 
licensed by governments [1|. Regulatory bodies have come to realize and announce that most of the time, large 
portions of certain licensed frequency bands remain unused. The intuitive intention behind secondary spectrum 
licensing is to efficiently increase the spectral usage of the network, while, depending on the type of licensing, 
not perturbing the primary users (higher priority users). Cognitive radio (CR) systems are seen as a candidate 
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prime solution that can significantly mitigate the current low spectral efficiency in the electromagnetic spectrum. 
A CR is defined as an intelligent wireless communication system that is fully aware of its environment and 
uses methodologies of learning and reasoning in order to dynamically adapt its transmission parameters (e.g., 
operating spectrum, modulation, and transmission power) to access portions of spectrum by exploiting the existence 
of spectrum holes left unused by a primary system. 

Cooperative diversity is a recently emerged technique for wireless communications that has gained wide attention 
(2). Most of the work on cooperative communications has concentrated on the physical layer (PHY) aspects of 
the problem by improving transmission parameters at the PHY. In a wireless communication network with many 
source-receiver pairs, cooperative transmission by relay nodes has the potential to improve the overall network 
performance. 

In J5), the authors proposed the relaying scheme shown in Fig. Q] where the secondary transmitter is used as a 
relay for the undelivered packets of the higher priority user (primary user). Also, they suggested the use of admitting 
parameter to control the relaying fraction. In [4|, an extension of the problem with multiple secondary transmitters 
acting as relays for the undelivered packets of the primary user was proposed, with and without opportunistic 
sensing scheme. In addition the authors of (4) considered priority in transmission is given to the relaying queues. 
In [5 1, the authors assumed that the cognitive user transmitter will be allowed to use the channel if the primary user 
not using the spectrum, and that the priority of transmission is given to the relaying queue. Also, it was assumed 
that the secondary will decide to relay a certain fraction / (adaptive admission control parameter) of the undelivered 
packets of the primary user to minimize the secondary user delay subject to a power budget for the relayed primary 
packets. In (6], the authors characterized fundamental issues in a shared channel where the users have different 
priority levels. Also, they characterized the stable-throughput region in a two user cognitive shared channel where 
the primary user has unconditional access to the channel while the secondary user transmits its packets with some 
access probability p. In Q, the authors proposed a cluster of secondary users helping the primary user with a single 
relaying queue accessible by all the secondary users. In fl8], the authors considered a network with two primary 
users and one secondary user relays their undelivered packets in the free time slots. In [9|, a multiple primary users 
and one secondary user capable of relaying was considered. 

We consider a slotted time-division multiple-access (TDMA) framework in which each time slot is assigned only 
to one terminal, i.e., orthogonal multiple access. If the primary system does not have a packet to transmit in its 
time slot, then this time slot is not utilized. These unutilized time slots are wasted channel resources that could be 
used to enhance the system performance and spectral efficiency. The cognitive system will be allowed to utilize the 
unused channel resources. 

In this work, we propose a novel cognitive system with the concept of cooperation. In the cooperative secondary 
transmitter-receiver system (CSTR), the cognitive system with its transmitter receiver pair tries to utilize the periods 
of silence of the primary transmitter in order to increase the reliability of communications against random channel 
fades and to allow the secondary user to utilize the channel effectively. In particular, when the secondary system 
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senses the channel for empty time slots, the slots are then used to help the primary system, and/or to allow the 
secondary user's packets to be served. In the non-cooperative secondary receiver with no priority queueing (NCNPQ) 
system, only the secondary transmitter helps the primary user to deliver its undelivered packets by certain access 
probabilities assigned to each queue. In the non-cooperative secondary receiver with priority queueing (NCPQ) 
system, the traffic of the secondary transmitter will not be served until the relaying queue empties. 
We make the following contributions in this paper. 

« We propose a new cooperative system (CSTR), which to the best of our knowledge, has not been proposed 
before in a cognitive networks where the receiver of the cognitive user helps the secondary transmitter to 
utilize the spectrum. 

< The main differences between CSTR system and [4| are as follows. In |4|, the authors considered multiple 
users each of them has its own arrival rate and a priority in transmission was given to the relaying queue. 
Also, the authors provided only an inner bound on the secondary network stability region. In this work, we 
investigate the stability region of a secondary system in presence of a primary system. We propose queues' 
access probabilities, so there is no priority in transmission given to any of the queues. Also, we consider 
another optimization parameter which is the priority of keeping (keeping factor) the undecodable primary 
packet at the primary receiver. We provided inner bound and outer bounds on the stability of the primary- 
secondary network, i.e., the maximum stable-throughput of the secondary transmitter given certain average 
arrival rate to the primary transmitter. 

« We characterize the optimal resources (e.g., adaptive acceptance factors, and access probabilities assigned to 
each queue) used by the secondary transmitter and its receiver, which maximize the stable-throughput region 
of the network. 

« We study the stability of two known systems and we characterize the optimal resources which maximize the 
maximum stable-throughput region for each system. 
Our theoretical reveal significant performance gains of the CSTR and NCNPQ systems over NCPQ system, and 
the significant improvement in performance of the CSTR system over NCNPQ system. 

The rest of the paper is organized as follows. In the next section, we describe the system model considered in 
this paper. The stable-throughput regions of the proposed systems are considered in Section [[II] In Section HVI we 
provide some numerical results, and finally, conclusions are drawn in Section [V] 

II. SYSTEM MODEL 

We consider the cognitive relaying system depicted in Fig. [2] where a secondary (cognitive) user is utilizing 
the spectrum resource whenever the primary user is idle. We assume that the secondary transmitter will be able 
to sense the channel every time slot to check whether the primary user is idle or not. The cognitive system will 
be able to send a packet each time slot during the idle sessions of the primary user. The main assumptions of the 
system model at both the MAC and PHY layers are given in this section. 
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A. PHY Layer Assumptions 

Denote the primary source as ps, the primary receiver as pd, the secondary source (transmitter) as ss, and 
the secondary receiver as sd. Let h* k denotes the channel gain between node j and node k at instant t, where 
j,k £ {ss, sd,ps,pd}, and it is distributed according to a zero mean circularly symmetric complex Gaussian 
random variable with variance cr? fc , i.e., CJ\f(0, cr? fc ). Channel gains are independent and each link is perturbed by 
complex additive white Gaussian noise (AWGN) zj k with zero mean and with variance N Q and independent for 
all links. In this paper, we consider multipacket reception (MPR) channel model which can capture the effect of 
interference, and fading at the PHY layer. Packets could survive the interference caused by concurrent transmissions 
if the received signal-to-interference and noise ratio (SINR) exceeds the threshold required for successful decoding 
at the receiver, i.e., for the link between node j and node k the probability of successful reception is given by 
P c j,k = Pr{SINR > jth, k}> where the superscript c denotes concurrent transmission event (see Appendix A). In 
case of no concurrent transmissions, successful transmission is achieved when the instantaneous signal to noise 
ratio (SNR) at the receiving node, employing the transmission power Pj, given by k \ 2 j^-, exceeds a predefined 
threshold ~fth,k, i- e -> the receiver is assumed to be able to decode the received message with negligible probability 
of error. The signal received at a receiving node k from a transmitting node j at time t can be modeled as 

yj,k = ^j,k X j + z ],k (1) 

where denotes the transmitted packet with average power Pj at time t. Given the channel model described above, 
the outage probability between node j and node k can be calculated as follows: 

Pr{O i , k } = Pj, k = PvUhl^Pj < N ^ th , k \ 

-~<th,k N o 

= 1 - e °'U p i (2) 

where Oj f. is the event that the link between node j and node k is in outage. The SNR threshold jth.k is a 
function of different factors in the communication system; it is a function of the application, the modulation, the 
signal processing applied at encoder/decoder sides, error-correction codes, and many other parameters [2 j. For more 
details regarding MPR the reader is referred to [6], [101- 11131 . 

B. MAC Layer Assumptions 

We assume that for the three systems the primary transmitter has a buffer Q p to store the incoming traffic 
packets, while the secondary transmitter has two buffers, Q s to store its own arrived traffic packets and Q ps to 
store the fraction of the undelivered packets of the primary user. We assume all buffers are of infinite length. 
We consider time-slotted transmission where all packets have the same size and one time slot is sufficient for the 
transmission of a single packet. The arrival processes of the primary and the secondary transmitter are assumed to 
be independent Bernoulli processes with mean arrival rates \ p and X s packet per time slot, respectively. Next we 
discuss the difference between the three systems. 
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1) Cooperative Secondary Transmitter-Receiver System (CSTR); In this system the secondary 
transmitter will accept a fraction of the undelivered packets of the primary user, i.e., it will decide to accept a 
controllable fraction to be admitted to its relaying queue in case of outage on the primary channel using an adaptive 
admission control parameter f s , and the secondary receiver has a buffer Q s d to store undelivered packets of the 
primary user in order to accept a fraction f s d of the undelivered packets of the primary user. We assume that for 
successfully decoded packets by the secondary transmitter and the secondary receiver the priority of keeping the 
packet is one of the optimization parameters of the system denoted by a binary value V, i.e., V = 1 if the priority 
is assigned to the secondary receiver, and V = if the priority is assigned to the secondary transmitter. We assume 
that the ACKs and NACKs messages initiated by the node with higher priority of keeping are sent earlier than the 
messages that are sent by the lower priority one, i.e., the node with priority of keeping transmits ACKs and NACKs 
from ti to T2 within the time slot and the other node transmits from T2 to T3. Note that the primary receiver sends 
the feedback messages over the period r Q to t\. The MAC layer is assumed to obey the following transmission 
scheme. 

• Assigns the priority of keeping the packet to the secondary transmitter or the secondary receiver at the 
beginning of the transmissions. 

• The primary transmits the packet at the head of its queue in its assigned time slot. If the queue is empty the 
time slot is free. 

• If a packet is received successfully by either the primary receiver, the secondary transmitter, or the secondary 
receiver the packet is removed from the primary transmitter's queue (the secondary transmitter or the secondary 
receiver needs to send an ACK if a packet is not decoded correctly by the primary receiver in this case). 

< If both the secondary transmitter and receiver decoded a packet correctly, the terminal which has the priority 
of keeping, will keep the packet and the packet will be removed from other node's buffer. 

• If a packet is not received successfully by the primary receiver, the secondary transmitter, and the secondary 
receiver, the primary transmitter retransmits this packet in its next assigned time slot. 

< At each sensed empty time slot, the secondary transmitter attempts to transmit a packet from its own queue 
with probability p s , it retransmits a packet from the relaying queue with some probability p sp , or it remains 
idel with probability pi and the secondary receiver attempts to retransmit the undelivered packets of the 
primary user with probability p s d, or to remain idle with probability 1 — p s d- 

• Packets could survive the interference caused by concurrent transmissions between the secondary transmitter 
and the secondary receiver, if the SINR exceeds the threshold required for successful decoding at the receiver. 

We assume that the overhead for transmitting the ACK and NACK messages is very small compared to packet sizes. 
The second assumption we make is that the errors and delay in packet acknowledgement feedback is negligible, 
which is reasonable for short length ACK/NACK packets as low rate codes can be employed in the feedback channel 
(2). In addition, nodes cannot transmit and receive at the same time, these transmission constraints are common in 
network systems where terminals are equipped with single transceivers J6). 
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2 ) Non-cooperative Secondary Receiver with Priority (NCPQ) and with no Priority (NCPNQ): 
In NCPQ and NCNPQ, the secondary transmitter is the only terminal that helps the primary user, it decides to 
accept a controllable fraction f s of the undelivered packets of the primary user. The difference between NCPQ and 
NCPNQ is that during periods of silence of the primary user, in the NCPQ, a priority in transmission is given to 
the relaying queue, while in NCNPQ some access probabilities are assigned to the secondary transmitter's queues, 
i.e., secondary transmitter transmits a packet from queue Q s with probability p s and from the relaying queue Q ps 
with probability p sp . The MAC layer is assumed to operate according to the following operating rules. 

• The primary transmits the packet at the head of its queue in its assigned time slot. If the queue is empty the 
time slot is free. 

• If a packet is received successfully by either the primary receiver or the secondary transmitter, the packet is 
removed from the primary user's queue (the secondary transmitter needs to send an ACK if one is not heard 
by the primary receiver in this case). 

• If a packet is not received successfully by both the primary receiver and the secondary transmitter, the primary 
transmitter retransmits this packet in its next assigned time slot. 

• At each sensed empty time slot, the secondary transmitter retransmits a packet at the head of its queue either 
from its own queue or from the relaying queue. 

III. Stability analysis of the systems 

A. CSTR 

Let us denote the queue sizes of the transmitting terminals at any time instant t by Q\, where i reads p for the 
primary transmitter's queue, s for secondary transmitter's queue, ps for the secondary transmitter's relaying queue, 
and sd for the secondary receiver's relaying queue. Then, Q\ evolves according to 

Ql+^iQl-U^ + A\ (3) 

where U\ is the number of departures in time slot t. A\ denotes the number of arrivals in time slot t and is a 
stationary process by assumption with finite mean £{»4*} = A^. The function (.) + is defined as (x) + = max(x, 0). 
We assume that departures occur before arrivals, and the queue size is measured at the beginning of the time slot 

an. 

A fundamental performance measure of a communication network is the stability of its queues. We are interested 
in the queues size. More rigourously, stability can be defined as follows |2), fl4l . 
Definition: Queue i S {p, s,ps, sd} is stable, if 

lim Pr{Q^ < y} = F(y) and lim F(y) = 1 (4) 

t— >oc y— >oo 

If the arrival and service processes are strictly stationary, then we can apply Loynes's theorem to check for stability 
conditions Q, 03). This theorem states that if the arrival process and the service process of a queue are strictly 
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stationary processes, and the average service rate is greater than the average arrival rate of the queue, then the 
queue is stable, otherwise the queue is unstable. 

The service and arrival rates of the nodes are as follows. For terminal i = p, given that the priority factor V = 1, 
i.e., the priority of keeping the packet is assigned to the secondary receiver, a packet can be served if either one 
of the four events is true: 1) The primary channel is in outage, the secondary receiver decides to accept the packet 
(with probability / st j), and the channel h ps sd is not in outage; 2) the primary channel is in outage, the secondary 
transmitter decides to accept the packet (which happens with probability f s ) and the secondary receiver decides not 
to accept the packet (which happens with probability 1 — f s( j), and the associated link h ps ss is not in outage; 3) 
the primary channel is in outage, the secondary transmitter and receiver both of them decide to accept the primary 
packet and both of them decode it correctly ^ or 4) if the channel between the primary source and primary receiver 
is not in outage, i.e., O ps pd is truj^. Due to the stationarity assumption of the channel gains, and using the outage 
probability formula ©, the probability of the event O pspd is given by Pr{0* s pd } = P ps , P d- The service process 
can be modeled as 



z4 = £i[4] 



(5) 



fc=i 



where 1[.] denotes the indicator function, and Ak, for k = {1,2,3,4}, are the events described above. From the 
above argument, it is clear that U p is stationary process and has a finite mean: 



fsdPps,sd (1 f sdPps,sd) f sP ^ 



ps.ss 



^{^p} — f^P — Pps,pd ~\~ Pps,pd 

where £{.} is the expected value. According to Loynes, stability of the queue i is given by: 

Hi > Aj. 

Therefore the condition on the stability of queue p is given by: 



Ap *C Pps.pd ~i~ Pps,pd 



fsdPps.sd ~t~ (1 JsdPps,sd)fsPp 



(6) 



(7) 



(8) 



If we take the priority of keeping factor into account, the general formula of the average service rate of the primary 
transmitter is given by: 

f^p Pps,pd Pps,pd *P ( fsdPps,sd ~\~ (1 f sdP ps,sd) f sP ps, 



ps,pd 



P 



ps : pd 



P\ fsPps,ss ~\~ (1 fsPps,ss)fsdPp 



fsd-Pps,sd 4" fsPps,ss fs-Pps,ssfsdPps,sd 



(9) 



'The primary packet will be buffered to the secondary receiver queue and dropped from the secondary transmitter queue due 
to the priority of keeping assigned to the secondary receiver. 
2 (.) denotes the complement of the event, i.e., x = 1 — x. 



s 



It should be mentioned that jj, p is independent of V. For terminal i = s, the service rate process can be modeled as 



{Ql = 0} p| A\ f) & ssM f]{Ql d - 0} + 1 {Ql = 0} f| Al f| d ssM [}{Q\ d ± 0} f) 



(10) 



where {Q* = 0} is the event that the primary queue is empty in time slot t, A f s denotes the event that in time 
slot t, the secondary transmitter assigned the channel to the relaying queue, i.e., the secondary transmitter will 
access the channel by the relaying queue (which happens with probability p s ), the secondary receiver is idle (which 
happens with probability p sd ), 0* s sd denotes the complement of the outage event between secondary transmitter 
and secondary receiver, {Q\ d — 0} is the event that the secondary receiver queue is empty, and A sd is the event 
that the secondary receiver is idle. The probability that the primary terminal queue is empty is given by: 



Pr{Q< = 0} = 1 - 



(ii) 



From the above argument, and the expression given in (fTTj) . it is clear that Ul is stationary process and has a finite 
mean: 



£{14} = Ms 
= (1 



sdPs 



Pr{Q* d = 0}+P sd Pr{Q* sd7 ^0} 



(12) 



According to Loynes the condition on the stability of the queue s is given by: 



X s < (1 - ^)P S 



sdPs 



Pr{% = 0}+P sd Pr{Q* sd ^0} 



(13) 



Consider now the relaying queue of the secondary transmitter Q ps . Given that the queue of the primary user is 
empty in a time slot t, a packet from queue ps can be served in either one of the three following events: 1) If the 
secondary receiver is idle (which happens with probability p sd ), and the channel between the secondary transmitter 
and the primary receiver is not in outage (9* s pd ; 2) if the secondary receiver accesses the channel (which happens 
with probability p sd ), and its queue is not empty, i.e., Q\ d ^ 0, and the complement of the event outage between 
the secondary transmitter and the primary receiver 0* S)P(2 is true; or 3) if the queue sd in time slot t is not empty, 
the secondary receiver accesses the channel (which happens with probability p sd ), and the complement of the 
event outage between the secondary transmitter and the primary receiver 0* s pd \T s d,pd- Mathematically, this can be 
modeled as: 

3 

Z& = X>[i3£] (14) 
fe=i 

where Ej,, for k = {1, 2, 3}, are the events described above. The expected value of the service process of the queue 
ps is given by: 



£{UL} = Mps 



MQld = 0} +P sd Pr{Qld ¥> 0} )Pss, P d +p sd Pr{Qld ^ o}p- 



~ss,pd 



\ (15) 

(1 - ^)Psp 
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Using Loynes conditions for stability and the expression of the outage probability in $Z§, the stability of the queue 
ps is given by: 

A. 



A ps < (1 



-)Psp 



(16) 



Consider now the secondary receiver queue Q sd . Given that the queue of the primary user is empty in time slot t, a 
packet from queue sd can be served if in a time slot t in either one of the six following events: 1) If the secondary 
receiver decides to access the channel (which occur with probability p s d), the secondary transmitter has no packets 
in any of its queues, i.e., Q\ = and Q* s = 0, and the complement of the event outage between the secondary 
receiver and the primary receiver O sd pd ; 2) if the secondary receiver decides to access the channel, the event that 
the secondary transmitter's queues are Q ps ^ and Ql — and it does not access the channel (which happens with 
probability p sp ), and the complement of the event outage between the secondary receiver and the primary receiver 
0* sd pd ; 3) if the secondary receiver decides to access the channel, the event that the secondary transmitter's queues 
are Q* s = and Q* ^ and it is not accessing the channel (which happens with probability p s ), and the complement 
of the event outage between the secondary receiver and the primary receiver 0* rf pd ; 4) if the secondary receiver 
decides to access the channel, the event that the secondary transmitter's queues are Q ps =^ and Ql ^ and 
it does not access the channel (which happens with probability pi = 1 — p sp — p s ), and the complement of the 
outage event between the secondary receiver and the primary receiver o\ d pd ; 5) if the secondary receiver decides 
to access the channel, the event that the queue Q* s ^ and the secondary transmitter accesses the channel (with 
probability p sp ), and the complement of the outage event between the secondary receiver and the primary receiver 
given a transmission between the secondary transmitter and the primary receiver 0* d pd \T ss ,pd\ or 6) if the secondary 
receiver decides to access the channel, the event that the queue Ql =^ and the secondary transmitter accesses 
the channel (with probability p s ), and the complement of the outage event between the secondary receiver and the 
primary receiver given a transmission between the secondary transmitter and the secondary receiver O sd pd \T ss ,sd- 
This can be modeled as: 

6 

i4d = ^[H] a?) 



fc=i 



where F{, for k = {1,2, 
queue sd is given by: 



. , 6}, are the events described above. The expected value of the service process of the 



£{Ut d } = fi sd 
= (1- 



— )Psd 



Pr{Q^ s = 0, Ql = 0} + p sp P,{Q ps ± 0, Q* = 0} + p s Pr{Q* s = 0, Q\ ^ 0} 



ftPr{Q* ^ 0,Qt ± 0} )P sdiPd + (pspMQls + 0} +p s Pr{Ql 0})P c sd: 



The condition on the stability of the queue sd is given by: 



pd 



(18) 



(19) 
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The arrival process to the queue ps can be described as follows. Given that V = 1, the primary user's queue 
is not empty, i.e., {Q P > 0}, the associated channel between the primary terminal-primary receiver is in outage, 
the secondary transmitter decides to accept the packet, and the channel between the primary transmitter and the 
secondary transmitter is not in outage, hence, the arrival to Q ps is either one of the following events: 1) The event 
that the secondary receiver decides to accept the packet from the primary user, and the associated channel between 
the primary transmitter and the secondary receiver is in outage; or 2) if the secondary receiver decides not to accept 
the packet. The process is modeled as 

2 

fc=l 

where Wt's are the events described above, and Pr{Q* > 0} = — . The process is stationary and the expected 
value of the arrival process is expressed as 

^■ps — Pps,pd{^ JsdPps,sd)fsPps,ss- (21) 
fJ-p 

Adding the priority factor, the mean arrival rate of the queue Q ps is given by: 

^ps — Pps,pd ( 1 P fsdPps,sd J fsPps,ss- (22) 

Mp \ / 

The arrival process to the queue Q s d can be described as follows. The event that the primary has packets, i.e., 
{Q p > 0}, the secondary receiver decides to accept a packet from the primary user W\ d , the channel between the 
primary user and the secondary receiver is connected, and the associated channel to the primary terminal-primary 
receiver is in outage. The process can be modeled as 



wu p|{Q* > o} n o ps , pd pi o psM 



(23) 



The process is stationary and the expected value of the arrival process to the queue Q s d is expressed as 



^{"^sdl — ^sd — fsdPps,pdPps,sd • (24) 

A*p 



If we involve V, the mean arrival rate of the secondary receiver queue is given by: 

^sd — Pps.pd f 1 *P ' JsPps,ss J fsdPps,sd- (25) 
V / 

Since the mean service rates at nodes s, ps and sd depend on each others queue size, these queues are called 
interacting queues, and consequently the rates of the individual departure processes cannot be computed directly. 
In order to overcome this problem, we utilize the idea of stochastic dominance J2), O, fl6l . ifTTl . which has been 
applied before to analyze interacting queues. Next, we discuss the stochastic dominance approach in order to obtain 
an inner bound on the stability region, then we construct two dominant systems, the inner bound on the stability of 
CSTR system is the union over the two dominant systems, and finally, we construct an outer bound on the stability 
region. 
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B. CSTR: Inner Bound 



To analyze the stability of the systems queues we resort to a stochastic dominance approach. Analyzing the 
stability of interacting queues is a difficult problem that has been addressed for ALOHA systems initially in [16|- 
fi"8l . Also, it should be mentioned that characterizing the stable throughput region for interacting queues with 
terminals greater than three is still a complex open problem. The basic idea of the stochastic dominant system 
technique is to decouple the interaction between queues and thus renders the analysis possible and easier. Also, 
it permits the characterization of the stability region as union of the stability regions of two dominant systems. 
It is done by forcing a node to access the channel in a consistent manner even when its packet queue is empty 
by letting the node to transmit dummy packets. Specifically, transmitting dummy packets is only aimed to cause 
invariant (constant) interference to the other node regardless of the emptiness of the queue and not counted as an 
actual packet transmission, i.e., the dominant systems and the original system are identical, except that, the packets 
successfully transmitted are not removed from the terminal's queue. 

In the stochastic dominance approach, we construct an appropriate hypothetical dominant system, which is a 
modification of the original system, that ensures that the queue sizes in the dominant system are, at all times, at least 
as large as those of the original system provided the queues start with identical initial conditions. Thus, the stability 
region of the new system inner bounds that of the original system. Furthermore, in the hypothetical system the 
queues are decoupled and not interacting anymore, thereby permitting the characterization of the stability region. 

1 ) First Dominant System : Let us denote by S the CSTR system. In the first dominant system designated 
as Si, where Q s d and Q s send dummy packets when their queues empty, and the Q. ps behaves exactly as it would 
in the original system S. Now, we can write down the service and arrival rates of the interacting queues, i.e., Q s , 
Q ps and Q s d as follows 

t's 



An 



(1 - - Pss,sd)PsPsd 

fj'p 



(26) 



ftps 



MP 



Psd^ 'ss.pd ~t~ PsdP C ss,pd 



The probabilities Pr{Q s = 0} and Pr{Q ps = 0} are given by: 



(27) 



Pr{Q s = 0} = 1- — 

Ms 

Pr{Q ps =0} = 

ftps 



ftsd z 



-(l-^)Psd 

MP 



ftps ftps / 



ftps 



ps : sd 



(28) 
(29) 

(30) 



According to the construction of the dominant system Si, it is easy to see that the queues of the dominant system 
are never less than those of the original system, provided they are both initialized identically (with the same initial 
conditions for queue sizes in both the original and dominant system). This is because, in the dominant system Si, 
the secondary receiver transmits dummy packets even if it does not have any packets in its queue, and therefore 
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interferes with secondary transmitter in all cases that it would in the original system. Therefore, given stability 
condition on the arrival rate for the queue Q s d if for some (A P ,A S ), the queues at all nodes are stable in the 
dominant system, then the corresponding queues in the original system must be stable. The portion of the stable- 
throughput region 1Z(S\), which is based on S\ is given by the closure of the rate pairs (A p , A s ) constrained by 
equations shown above as f s , f s d, p s , Psp, and p s d vary over [0, 1], and V varies over {0, 1}. This is essentially 
similar to the argument presented in J6], ifTTI . and is applicable here as well. To obtain the closure of the rate pair, 
we utilize the constrained optimization technique. For a fixed X p , the maximum stable arrival rate to the secondary 
transmitter's queue is given by the following optimization problem (as in Q, |6l): 



max 

Ps,Psp,Ppd,fsJsd,'P 

s.t. 



X s — Ms 

< p s ,p S p,Psd, fa, fsd < 1) V € {0, 1}, Ps+Psp < 1 
Xp < Hp, Xps < Mps, X s d < Hsd- 



(31) 



2) Second Dominant System: The second dominant system is designated as S2, where the secondary 
transmitter is the one that send dummy packets when any or both its queues are empty, and the secondary receiver 
behaves exactly as it is in the original system S. The service rates of the system's queues are given by: 

A,, 



Hsd 



(1 - —)Psd 

Mp 



PiPps,sd ~t~ ijpsp ~t~ Ps)P C ps,sd 



The probability that queue Q s d is empty is given by: 

MQsd = 0} 



Xsd 
Hsd 



Hs 



(1 



)Ps 



sd Ps 



1 -Psd- 



X, 



Hps = (1 



A p , 



[l -p s d — )Pss,pd +Psd'-^-i'' 

Hsd 



Hsd 

X s d 
1 

Hsd 



(32) 



(33) 



(34) 



Based on the construction of the dominant system S2, it is trivial to notice that the queues of the dominant system 
are never less than those of the original system, provided they are both initialized identically. This is because, in 
the dominant system 1S2, the secondary transmitter s transmits dummy packets even if it does not have any packets 
in its queues (transmits a dummy packet from queue Q s with probability p s and from the relaying queue Q ps with 
probability p sp ), and therefore interferes with secondary receiver in all cases that it would in the original system. 
Therefore, given stability condition on the arrival rates of the queues (Q s d, Qps)> if f° r eacn mean arrival rates pair 
(A p ,A s ) the queues at all nodes are stable in the dominant system, then the corresponding queues in the original 
system must be stable. The second inner bound for the stable-throughput region of the CSTR system, 1Z(S%), which 
is based on the dominant system 1S2 can be obtained by formulating a constrained optimization problem similar 
to that discussed above for the CSTR system, where we fix X p and maximize A s as f s , f ps ,Ps,Ps P and p s d vary 
over [0, 1] and V varies over {0, 1}. The maximum stable-throughput of both dominant systems is depicted in Fig. 
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[3] We used the following set of outage probabilities P pStP d — 0.4, P SSyS d = 0.2, P ps . ss = 0.08, P ss ,pd = 0.08, 
Psd,pd = 0.02, and P ps . s d = 0.02. The optimization problems of the first and second dominant systems are solved 
numerically using MatLab optimization toolbox. Since the problems are nonconvex, the solver produces a locally 
optimum solution. To increase the likelihood of obtaining the global optimum, the program is run many times, say 
10000 times, with different initializations of the optimization variables. 

C. CSTR: Outer Bound 

Here we provide two outer bounds for CSTR. 

1 ) First Outer Bound: An outer bound for the CSTR can be obtained by using the following identities 
ED: 

MQld = 0} +PsdMQld ^ 0} < Pr{Q* d = 0} + Pr{QLi ^ 0} = 1. (35) 

And 

(pr{Ql d = 0}+p sd Pr{Q* d + 0}) P ss , P d + p sd Pr{Q* d 0}P^ ss , pd < P ss , P d- 
In addition, based on Bayes' theorem 

Pr{a, b} = Pr{a|6}Pr{6} < Pr{6} or Pr{a, b} = Pr{%}Pr{a} < Pr{a} (36) 
where a and b are two event. Thus, we can upper bound the following quantities: 



Pt{QI s = 0, Ql = 0} < Pr{Q ps = 0}, Pr{Q* s ^ 0, Ql = 0} < Pr{Q£ s + 0}, 

Pr{Q* s = 0, Ql ± 0} < Pr{Q* s = 0}, Pr{Q* s ^ 0, Q\± 0} < Pr{Q£ s ^ 0} 
The mean service rates of the queues can be upper bound as follows 



(37) 



li. < (1 - — )^ ss , s dP s (38) 
Mp 

Mps < (1 - —)Ps P P ss , P d (39) 
A*p 



Therefore, the service rate of the queue Q s d is given by: 



Hsd < (l-—)p s 

MP 



((l-^) + P sp ^ + P s (l-^)+ Pi ^)p ps>sd+ ( Psp ^l +P s )p~o psM 

\ /^pS /^ps /^pS / \ ftps J 



(40) 



When the inequalities hold, the queues are not interacting anymore, therefore, we can obtain the outer bound by 
solving a constrained optimization problem to get the closure (A p , A s ). The optimization problem is given by: 

max A s = (1 - —)P ss ,sdPs- 

Ps,Psp,P P d,fs,f l! d,'P £tp 

S.t. < p s ,p spi p sd , fa, fad < 1, V € {0, 1}, Ps + Psp < 1- 

Ap < /Up, A ps < [l ps ,\sd < fisd- (41) 



14 



Denote the first outer bound as s[°\ the outer bound can be specified by the rates pair 



H(s[ o) ) =i (X P ,X S ) : A s < (1- 



Xr. 



Pps,pd ~t~ Pps,pd 



■)Pss,adP\ 



(42) 



fsd-Pp8,sd + (1 f s dPps,sd)ftPps,; 

where p*, p* sp , p* d , /*, /* d and are obtained by solving fiTY 

2) Second Outer Bound: Another outer bound which can be addressed analytically is obtained as follows. 
Using ©: 



jlp — Pps,pd ~l~ Pps,pd 



fsdPps,sd ~h fs-Pps,ss fsPps,ssfsdPps,sd 



^ Pps,pd ~\~ Pps,pd rn&X s Pps,sdi P 



ps,ss 



I*.<(1- —)P SS , S dPs < (1 - ^)P ss ,sd 

Mp Mp 



when the inequality d43l holds and applying Loynes's law 
A s < n„ < (1 



Pvs.vd ~"H Pvs.vd max ^ Pps : sdi Pps,ss 



)Ps3, 



sd 



A, < (1- 



Pps,pd H~ Pps : pd rXlSX s Pps^sdi P 



ps,ss 



^)P ss,sd 



(43) 
(44) 

(45) 
(46) 



Denote the second outer bound as (S^ , the outer bound can be specified by the rates pair 



/^(^ o) H(A p ,A s ) A 



X v 



p 



ss.sd 



p. 



ps,pd ~r ±ps,pd 



Pnti.nd. niRX S P ps^sdi Pps,ss 



< l 



(47) 



The outer bound, <S (o) , of the CSTR is the intersection of the two outer bounds, i.e., 7e(5 (o) ) = TZ(S[ o) ) f| K(S. 



4 0) ) 



Note that K(s[ o) ) is contained inside K{S y 2 °>), i.e., 1Z{S\°>) f| Tl(S 2 0> ) = TZ{S\°'). The inner and the outer bounds 
of the CSTR are depicted in Fig. g] 



;(°) 



(oh 



D. NCNPQ 

In this subsection, we investigate the stability region of a cooperative cognitive user with no priority assigned 
to the relaying queue and with adaptive acceptance factor. The secondary transmitter will use the idle time slots of 
the primary user to transmit a packet from the relaying queue with probability p sp or to transmit a packet from its 
own traffic with probability p a = 1 — p sp . The arrival and service rates of each node are given by: 

(48) 
(49) 



X 



f sPpS.pdP pS,S8 



ftp Pps.pd ~\~ fsPps,pdPps,s 

ftp 

f^ps — Psp(l ^)Pss,pd 
A*p 

Us = Ps(l - — )Pss,sd- 



(50) 
(51) 
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The closure of the pair of arrival rates (A p , A s ), where we fix X p and maximize A s as f s and p s vary over [0, 1]. The 
stable-throughput region of the system can be obtained by formulating a constrained optimization problem similar 
to that discussed above for the CSTR system. The optimization problem is expressed as: 

max A s = p s (l - — )P ss ,sd 

fs ,Pa flp 

s.t. Ap ^ fip 

0< fs < 1,0 <Ps < 1 

Aps 5= Mps (52) 

Arranging the above optimization problem as follows: 

Ps{P ps ,pd + f s K - Ap) 
max ^=^- — , 

f"P> Pps,pd + fsK 

S-t. —J S K ^ Pps,pd hp 

o< /. < 1,0 <p s < 1 



fs P SS.pd , P SS.pd srr*\ 

Ps T /s ^ t7\ (53) 



P ps ,pd + f s K - Ap KX P ~ K\ p 
where K = P ps , P dPps,ss- The optimization problem is non-convex. It can be solved via a one-dimensional grid 
search over the optimal value of f s . Fixing f s , we have the following optimization problem: 

max p s 

Ps 

S.t. f s > X P-^ d 

K 

0<fs<l,0< Ps <l 

Ps<l-= f - — J^- (54) 

Pps.pd ~\~ fsK Ap Pss,pd 

The optimization problem is convex and can be readily solved. The optimal value of p s is given by: 

* i f s kx p 

P s = 1 — = = — (55) 

Pps,pd ~t~ fsK Ap P ss,pd 

* f s K\ p 
p* = = = — — (56) 

pps.pd i j s-** ^p l ss.pd 

The stability of the NCNPQ system is given by: 

ft(NCNPQ) =((A p , X s ) :X s <p*Jl- = — ^ = )T' A. (57) 

I \ *ps,pd \ J s f ps,pd-Lps,S8 / J 

The optimal values are drawn in Fig. 
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E. NCPQ 

In this subsection, we investigate the stability region of NCPQ. The service and arrival rates of each queue in 
the system is given by 151 : 

f^p — Pps.pd ~t~ fsPps.pdPps.ss (58) 

Aps fsPps.pdPps.ss (59) 

lip 

lips = (l-^)Pss,pd (60) 
MP 

Us = (1 - ^)(1 - ^)Pss M - (61) 

flp flps 

Stable-throughput region of the system can be obtained by formulating a constrained optimization problem similar 
to that discussed above for the CSTR system, where we fix A p and maximize A s as f s varies over [0, 1]. 



max A s = (1 - ^)(1 - ^)P ss>sd 

Js \Xp ftps 

s.t. Xp (.ip 

< f s < 1 

A ps < lips (62) 
After some mathematical manipulation the optimization problem reduces to: 



mm 



l + fsj^ 



/• 1 + fs 



K 



s.t. u 7 '" 



A" 

< / s < 1 

/s(Ap — P ss ,pd) 5: (Pps.pd ~ Ap) , P (63) 

The constraints in the above optimization problem are linear. However, the objective function is not convex. 
Nevertheless, it is linear fractional, and hence, quasi-convex (quasilinear) |20|. The bisection algorithm can be 
employed to solve this kind of optimization problems iteratively by posing it as a sequence of feasibility problems 
where each problem is a linear program. In our case, it can be readily solved as in Appendix B. It can be noted 
that if Pss.pd = Pps,pd the above optimization problem will be reduced to a feasibility problem ||20l . 

The optimal value of the acceptance factor f s according to the channel quality between of the primary link and 
primary-secondary link is given by: 

• If Pps.pd ^ Pss.pd 

' [0,1] if A p = 0; 

1 if < A p < P ss ,pd • 



f*s 
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• ^ Pps,pd ■* 'ss.pd 

r= [ [0,1] ifA p = 0; 

[ if < Ap < P ps ^ pd . 

• If Pps 7 pd — P ss,pd 

The problem is reduced to a feasibility problem. However, we will obtain set of fractions f s that satisfy the 
domain of the objective function and the constraints. Thus the optimal set of f s is given by: 



A* = [0, 1] (64) 

Proof: See Appendix B. ■ 
Therefore, the secondary transmitter according to the channels quality will choose the optimal value of the admitting 
factor f s . 

The stability region of the NCPQ system is given by: 

^(NCPQ) = S(X P , A.) : A. < (1 - \ )(1 - =^-= f ^ _ )P sa ,A (65) 

k *ps,pd T" Js ss,pd * ps.pd ' Js P ) 

Proposition 1. 77ie stability region of NCPQ is contained inside that of CSTR and NCNPQ. In other words 

ft(NCPQ) C ft(NCNPQ) C ft(CSTR) (66) 

Proof: Since the NCPQ is achievable by the CSTR and NCNPQ system via controlling the optimization 
parameters (setting some of them to zero), therefore, its stability is achievable and contained inside both of them. 
On the other hand, NCNPQ is achievable by the CSTR via controlling the optimization parameters, hence its 
stability is contained inside CSTR. ■ 
Next, we provide some numerical results to reveal the above proposition. 

IV. Numerical Results 

In this section we provide the solution of the optimization problem considered in this paper. A comparison 
between the proposed systems is depicted in Figs. [6] and [8] for different values of outage events and probability 
of reception in case of concurrent transmissions. The corresponding optimal values of the optimization parameters 
for the CSTR are drawn in Figs. [7] and |9] respectively. It should be mentioned that if the link of the primary is 
always in outage with probability 1, i.e., there is no direct link between the primary transmitter and receiver (the 
primary user's packets will never be served if there is no cooperation, i.e., [i v — P ps ,pd — 0) and/or the link 
between the primary receiver and the secondary transmitter is in outage most of the time, the CSTR performance 
will significantly overcome the maximum stable throughput of the NCPQ and NCNPQ systems. This remark can be 
seen in Figs. [T0land[TT1 The optimal values for CSTR to achieve the boundary of Fig. QT|are depicted in Fig. [12] 
respectively. The mean arrival rate of the secondary transmitter is a monotonic decreasing function with the mean 
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arrival rate of the primary transmitter, A p . This is because as the the primary arrival rate increases the probability 
of the primary queue to be empty vanishes and the secondary queue will not be able to access the channel. 

V. Conclusion 

In this paper, we have addressed the impact of cooperative cognitive system at the multiple-access layer. We 
investigated the maximum stable-throughput of three systems, i.e., CSTR, NCPQ, and NCNPQ systems. In CSTR, 
the cognitive system (transmitter and receiver) senses the channel for idle channel resources and exploits them to 
relay the undelivered packets of the primary user or to use the channel resources to serve its own traffic. In NCNPQ, 
the secondary transmitter uses the silence periods of the primary user to serve a packet from one of its queues. 
In NCPQ, the secondary transmitter assigned higher priority to the relaying queue, if the relaying queue is empty, 
it transmits a packet from its own queue. The results reveal the significant performance of CSTR and NCNPQ 
systems over NCPQ. The maximum stable-throughput of the NCPQ is a subset of the maximum stable-throughput 
of NCNPQ, and the maximum stable-throughput of the NCNPQ is a subset of the maximum stable-throughput of 
the CSTR system. 

Appendix A 

The outage probability between terminal i and terminal j given that there is a transmission between terminal k 
and terminal j 0, is given by: 

^=Pr {0i ,|T fc , } = ^ pJmIn. ^ 



1 



1 



1 , Pklth,j a k,j 

where {Tk.j} is the event that the terminal k is transmitting a packet to terminal j. One can obtain a relationship 
between the two events {Oij\Tkj} and {Oij} as follows: 



p%3 = — /; j j2 • Pij (67) 

£V7thJ °k,j 
Pi 



1 -| -r- 



where the inequality of (e) holds to equality if Pf. = 0, i.e., no concurrent transmission. 

Appendix B 

In this Appendix, we provide the solution of the following optimization problem: 

1 + xPx 
mm — 

x 1 + xpi 

S.t. X > c 



< d (68) 
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where c,d,j3\ and /3 2 are constants. If /?i > (3 2 , then to minimize the objective function we need to choose the 
lowest value of x, i.e., x = c. If (3\ < #2, then to minimize the objective function we need to choose the highest 
value of x, i.e., x = d. On the other hand, (i\ — 02, the optimization problem reduces to: 

min 1 

X 

S.t. X > c 

x<d (69) 
which is a feasibility problem. The optimal set is x* = [c,d]. 
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. 1. Primary and secondary links of the NCPQ system and NCNPQ. 
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Fig. 3. Stability region of the dominant systems, the inner bound for the stability of the CSTR system is 1Z(Si) \J H(S2). 
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Fig. 5. The optimization parameters of the NCNPQ system: P ps . p d = 0.4, P SSjS( i = 0.2, P ps , S s = 0.08, P S s, P d = 0.08, 
P s d,pd = 0.02, and P ps , s(i = 0.02. 
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Fig. 6. Comparison between the maximum stable-throughput of the CSTR, NCPQ, NCNPQ, and the no cooperation case: 
P P s, P d = 0.4, P ss , sd = 0.2, P p3 ,ss = 0.08, P ss ,pd = 0.08, P sd , pd = 0.02, P ps , sd = 0.02, P c ss , pd = 0.2640, and P c sd , pd = 
0.2160. 
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Fig. 8. Comparison between the maximum stable-throughput of the CSTR, NCPQ, NCNPQ, and the no cooperation case: 
P P s, P d = 0.4, P ss , sd = 0.2, P ps , ss = 0.08, P ss , P d = 0.08, P sd , pd = 0.02, P ps , sd = 0.02, P c ss , P d = 0.54, and P c sd , pd = 0.51. 
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Fig. 9. The optimization parameters of the CSTR system: P pStP d = 0.4, P ss . s d = 0.2, P ps ,ss = 0.08, P S s, P d = 0.08, 
Psd, P d = 0.02, P ps , sd = 0.02, P c a3:Pd = 0.54, and P c sd , P d = 0.51. The optimal value of the priority factor P* = 1. 
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Fig. 12. The optimization parameters of the CSTR system: P pStP d = l,P S s,sd ~ 0.2, P pSjS3 = 0.8, P ss ,pd = 0.8, P s d 
0.02, P ps , sd = 0.02, P c as , vd = 0.54, and P% d , v d = 0.51. 



